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We study a new system in which electrons in two dimen- 
sions are confined by a non homogeneous magnetic field. The 
system consists of a heterostructure with on top of it a su- 
perconducting disk. We show that in this system electrons 
can be confined into a dot region. This magnetic anti-dot has 
the interesting property that the filling of the dot is a dis- 
crete function of the magnetic field. The circulating electron 
current inside and outside the anti-dot can be in opposite di- 
rection for certain bound states. And those states exhibit a 
diamagnetic to paramagnetic transition with increasing mag- 
netic field. The absorption spectrum consists of many peaks, 
some of which violate Kohn's theorem, and which is due to the 
coupling of the center of mass motion with the other degrees 
of freedom. 
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I. INTRODUCTION 

Quantum dots have been the subject of botk theo- 
retical and experimental research in recent years.B They 
have been successfully created experimentally by apply- 
ing lithographic and etching techniques to impose a lat- 
eral structure onto an otherwise two-dimensional elec- 
tron system. These structures introduce electrostatic 
potentials in the plane of the two-dimensional electron 
gas (2DEG), which confine the electrons to a dot region. 
The energy levels of electrons in such a quantum dot are 
fully quantized like in an atom, and therefore are also 
referred to as artificial atoms. In such electrically con- 
fined quantum dots the confinement potential can well 
be represented by a parabolic potential. 

In the present paper we study a new quantum dot sys- 
tem which is different from the usual quantum dot sys- 
tem: 1) the electrons are confined magnetically, 2) the 
confinement potential is inherently non-parabolic, and 3) 
the dot contains a finite number of electrons where the 
filling of the dot is a discrete function of the strength of 
the confinement (magnetic field). j-. 

This magnetic anti-dot can be realizeda by deposit- 
ing a superconducting disk on top of a 2DEG. When a 
homogeneous magnetic field is applied perpendicular to 
the 2DEG, the magnetic flux lines are expelled from the 
superconducting disk due to the Meissner effect, which 



results in an inhomogeneous magnetic field profile in the 
2DEG. Note that this problem is related to the problem 
of type II superconductors on a heterostructure, where 
flux lines penetrate the 2DEG.B Between the flux lines 
regions of zero magnetic field are present. Here we have 
the inverse situation in which we have a uniform magnetic 
field except for a local dot-like region where there is no 
magnetic field present. Bound states are now possible in 
the zero magnetic field region. 

Such a system was discussed by Peeters et al. in Ref. ^ 
where preliminary results were presented for the energy 
levels of such a magnetic anti-dot. Here we elaborate 
on this system and give a more detailed and complete 
study of the bound states of such a system and calculate 
also other properties of this system. In Ref. || the Hall 
and bend resistance resulting from such a system was 
discussed in the ballistic regime and in Ref. ^ the diffusive 
transport of such a magnetic anti-dot placed on top of a 
Hall bar was studied. 

Solimany et alu studied a limiting case of the present 
system in which the magnetic field was B{p) = Ba5{p—a) 
and they solved the classical and quantum mechanical 
equations for a majsnetically confined quantum dot, and 
recently Sim et ala investigated the formation of mag- 
netic edge states along with the corresponding classical 
trajectories. Although both groups referred to this the- 
oretical system as a magnetic dot, we think the name 
'magnetic anti-dot' is more appropriate in this case. Here 
we start from the experimental realizable system, i.e. a 
superconducting disk on top of a 2DEG, and calculate 
in Sec. II the non-homogeneous magnetic field profile in- 
duced in the 2DEG. We find that the magnetic field pro- 
file is different from the one assumed in Refs. and ||. 
In Sec. Ill the energy spectrum of electrons in the 2DEG 
near this magnetic field profile is calculated and com- 
pared with those of two circular model magnetic anti-dot 
systems. The filling of the magnetic anti-dot is calcu- 
lated as function of the strength of the confinement. In 
Sec. IV we discuss the electron probability current and 
the induced magnetic moment. The optical absorption 
spectrum, i.e. frequencies and oscillator strengths, are 
obtained in Sec. V. The results are summarized and our 
conclusions are presented in Sec. VI. 
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II. MAGNETIC FIELD PROFILE 

The system we have in mind is shown schematicahy in 
the inset of Fig. |l|. We have a high mobihty heterostruc- 
ture with on top of it a superconducting disk placed in 
a homogeneous apphed magnetic field B = {0,0, Ba). 
For convenience, we consider a very thin superconduct- 
ing disk which is perpendicular to the magnetic field and 
a distance z-above the 2DEG. Because of the symmetry 
of the system, we use cylindrical coordinates (p, ip, z) 
and measure all lengths in units of the disk radius a. We 
solve the following magnetostatic equations outside the 
disk 



V B 

V X B 



= 0, 
= 0, 



(1) 



with the conditions that = due to symmetry con- 
siderations, and i?2|p<i = on the disk surface (z=0). 

This problem is equivalent to the problem of liquid flow 
around a disk and was solve d in Ref. |9|. Ther e the oblate 
spheroidal coordinates p = \J (^^ + 1)(1 — rf'), ip and z — 
^rj were introduced, where 0<^<cxd,0<77<1 and 
< if < 2tt which have t he corresp onding scale factor s 

\ = + The surface ^ = is a disk 

of radius 1 in the xy-plane with center at the origin. 
According to that solution, the magnetic field under such 
an infinitesimal thin disk can be represented as B = — V$ 
where 



1 — ^ tan 



(2) 



Using cylindric symmetry and B = V x A, we arrive at 
the single non zero component of the vector potential 



— tan 



, (3) 



from which we obtain the perpendicular component of 
the magnetic field profile (or from Eq. (Q) using B — 



P 



2Bn 



(4) 



These results can be easily converted back to the 
cylindrical coordinate system [p, if, z) by inserting 
^2 ^ (l/2)[v/(r2 - 1)2 +4z2 + (r2 - 1)] and r/^ = 
(l/2)[v/(r2 - l)2 + 4z2 - (r2 - 1)] with r"^ = p^ + 7? in 
the above equations. 

The resulting magnetic field profile is shown in 
Fig. [^(a) for different values of the distance (z) between 
the 2DEG and the superconducting disk. Notice that 
the magnetic field under the disk is very small due to the 
Meissner effect, while far from the disk it becomes equal 
to the external magnetic field strength Ba- At the edge 
of the disk there is an overshoot of the magnetic field 



strength, which becomes larger with decreasing value of 
z. When the 2DEG is further away from the supercon- 
ducting disk, there is almost no overshoot and B^^ grad- 
ually increases away from the center. Notice that in the 
latter case the magnetic field is nonzero and consequently 
the model systems of Ref . |^ and Ref. || can not be real- 
ized by using a thin superconducting disk. 

Before discussing the electron states in such a mag- 
netic field profile, we will first consider two model sys- 
tems which correspond to two extreme situations, but 
which contain the essential physics of the problem. The 
two models we consider are defined by the following pro- 
files: 1) Bz{p) — BaO{p — 1) with corresponding vec- 
tor potential A^{p) = [Ba{p - l/p)/2]e{p - 1), and 
2) = {Bap/2)0{p — 1), which results into a mag- 
netic field profile with a delta function overshoot B^ = 
Ba6{p - 1) + {Bap/2)6{p - 1). The first model was used 
as the magnetic field profile in the magnetic anti-dot of 
the previous paperaJu. The second model contains a 
magnetic overshoot at the edge of the magnetic dot (see 
also Ref. ^), which models the real system for z/a <C 1. 
This overshoot is due to the de-magnetization effects of 
the superconducting disk. The vector potential profile 
of these two models is shown in Fig. |l|(b) together with 
the one resulting from the superconducting disk for two 
different values of the set back distance z. 



III. THE ENERGY SPECTRUM 

In order to calculate the wavefunctions and their 
corresponding energy, we have to insert the expres- 
sion for the spatial dependent vector potential (A) 
into the momentum operator p — » p + (e/c)A which 
results in the Schrodinger equation: — (l/2me)[?iA -|- 
(ie/c)A(r)]25'(r) = E'^{y). Due to the cyhndrical sym- 
metry of the problem, the wave function can be presented 



^{p,p>) = 



1 



27r 



(5) 



and the problem reduces to solving the ID radial equa- 
tion 



^ + iiL + 2[i?-U(p)]^i?(p) = 0, 
dp'^ p dp 



where the effective potential reads 



V{p) 



AM 



m 
P 



(6) 



(7) 



We solved this eigenvalue problem numerically, using the 
Newton iteration technique and subjecting the solution 
to the following boundary conditions: R{p 0) = yol™' 
and R{p oo) — 0. 

The numerical results for the energy spectrum are 
shown in Fig. ||(a) for the model without overshoot, in 
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Fig. ^(b) for the model with overshoot and in Fig ||(c) 
for the real profile in case of zja = 0.01. The different 
energy levels are labeled with the corresponding quan- 
tum numbers (n,m). We found it convenient to express 
the energy in units of i?o — Ti' jmeO^ and the applied 
magnetic field Ba in units of Bq = ch/ea^. These units 
are related to the problem of a particle in a box. For 
example for a = 0.01, 0.1, 1, 10/im we have respectively 
Eo = 7.63 X 10-1, 7.63 x lO'^, 7.63 x 10"^, 7.63 x 10"'^ 
meV and Bq = 6.6 x 10"*, 6.6 x 10^ 6.6, 0.066 Gauss. 
From Fig. ^ we notice that for small magnetic fields the 
energy is linear in Ba and in fact we recover the Lan- 
dau levels En,m = fiujc{n + {\m\ + m)/2 + 1/2) for an 
electron in a homogeneous magnetic field. The reason is 
that for small magnetic field we have Ib/o, ^ 1, where 
Ib = \/fic/eB is the magnetic length. So the electron 
wavefunction is spread out over a large region and most of 
its probability is in the p/a > I region where B{p) « Ba- 
in the opposite case of a large external magnetic field we 
have Ib/o, <ti 1 and the electron wavefunction is localized 
in the dot region where there is no magnetic field present. 
The problem is then similar to the one of an electron in 
a circular dot in the absence of a magnetic field and we 
recover the discrete energy levels which are determined 
by the zeros of the Bessel function: J|^|(fc) = and are 
plotted on the right hand side of the figure. In this case 
the energy scale is Eq — /meO^ like for a particle con- 
fined in a ID box. 

Although the limiting behavior of the spectrum for 
Ba — > and Ba — > oo is very similar in all three cases 
the intermediate behavior turns out to be very different. 
Indeed we see that the energy levels of the model 

without overshoot increase slower than linear with in- 
creasing Ba- This is not the case when there is a magnetic 
field overshoot where there are a number of energy levels 
which for small/intermediate fields have a superlinear be- 
havior as function of Ba- This distinct behavior is made 
more visible when we plot the energy in units of HlOc as 
function of 1/Ib ~ Vb2 as is done in Fig. ^. This dif- 
ferent behavior between the two cases can be understood 
by looking at Fig. H, where the radial part of the electron 
wavefunction (n,ni)=(0,-l) is plotted for various values 
of the magnetic field strength Ba in case of the model 
with magnetic overshoot. For m < the wavefunction 
exhibits a maximum at p — p* > 0. With increasing 
1/Ib this maximum shifts towards the center of the dot. 
In case there is an overshoot in the magnetic field pro- 
file and when the maximum of the electron wavefunction 
is situated near p/a = 1, the electron energy will be in- 
creased which results in the local maximum in the energy 
as shown in Fig. ^ . From this interpretation it is easy to 
understand that when m < the maximum in iJ^^m/TiWc 
shifts to larger 1/Ib with increasing \m\. 

This behavior of the electron energy has important 
consequences for the filling of the dot. Outside the quan- 
tum dot the magnetic field is Ba and the electron lowest 
energy state is huJc/'2. An electron will only be situated 



in the dot when its energy is lower than in the region 
outside the dot. From Fig. ^(a) we note that for a dot 
without magnetic field overshoot at its edges there are 
an infinite number of states, i.e. the states with m < 
for n = 0, which have an energy less than huJc/'2 and 
consequently the electrons will be attracted towards the 
dot. For the system with magnetic overshoot the situ- 
ation is totally different. The |0, 0) state has an energy 
below huJc/2 and two electrons (two because of spin) will 
be able to occupy the dot. When we try to add more 
electrons to the system we see that for small 1/Ib the 
electrons will prefer to sit far away from the dot region 
where they have a lower energy, i.e. fiUc/'i- Thus in 
this situation the electrons are repelled by the anti-dot. 
Increasing the magnetic field will bring the (0,-1) level 
below hojc and then two more electrons will be attracted 
towards the magnetic anti-dot. This discrete filling of 
the dot is shown in Fig. |^ for the model with magnetic 
overshoot (solid curve) and for the superconducting disk 
case (dashed curve) with z/a — 0.01. 

Including the real magnetic profile does not alter our 
conclusions qualitatively. This is shown for z/a = 0.01 
in Fig. ||(c) and Fig. ^(c). We obtain some kind of inter- 
mediate behavior between the two model systems. The 
discrete filling of the dot is still present, which can be 
inferred from Fig. ^(c). Nevertheless, the exact position 
at which the number of electrons jump to higher values 
is a function of the exact magnetic profile, and in partic- 
ular depends strongly on the sharpness of the magnetic 
overshoot. This is why z/a has to be very small, which 
can always be experimentally achieved, by making the 
superconducting disk large enough. 

IV. THE PROBABILITY CURRENT 
DISTRIBUTION AND THE MAGNETIC 
MOMENT 

The probability current distribution of the eigenstates 
is also different from the usual quantum dot case. When 
we rewrite the wavefunction as ^{v) — a(r)e*^^""\ the 
probability current is given byli3 

J(r) - (l/m,)a2 [nve(r) + (e/c)A(r)] , (8) 

where the first term is the well-known circular current 
and the second term is due to the magnetic field. For 
bound states the current vector has only an angular 
component J(r) — J^{r)e^ which is independent of ip. 
The probability current distribution Jip(p) and the cor- 
responding radial distribution function of bounded states 
with n = and different m-values (indicated in the fig- 
ure) are plotted in Fig. || for the three different cases. 
For TO = no current flows inside the dot except for the 
realistic magnetic field profile (Fig. ||), although it is very 
small. The to > states have a positive circular current 
which is larger with increasing overshoot of the magnetic 
field at the edge. For to < the current is circulating 
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in the opposite direction and the magnitude increases 
with increasing overshoot. The maximum of the cur- 
rent profile moves closer to the edge with increasing |m|. 
Sufficiently outside the dot region the circular current is 
positive, irrespective of the value of m. With the delta 
overshoot this is also true near the outside edge of the 
dot where the current distribution exhibits a discontinu- 
ous behavior. In the latter case the current intensity is 
a uniform decreasing function of the distance p > a. For 
the case without overshoot, or when we have a contin- 
uous overshoot, the current is negative near the outside 
edge for to < 0. These results can be understood from 
classical trajectories of magnetic edge states circulating 
clockwise (m < 0) or counterclockwise (to > 0) along 
the boundary region of the magnetic anti-dot without 
overshoot (Cfr. Ref. |^) and which are shown in Fig. ^ 
for to = —1, TO = and to = 1. For completion, we 
also included Fig. which is the same as Fig. ^, but for 
the states with to = 1 and different n-values (indicated 
on the figure). Now the radial part of the wavefunction 
has n-nodes which results in zeros in the circular current. 
Note also that the current becomes much more strongly 
peaked near p « 0. The amount of current outside the 
dot also increases with increasing n. Because to = 1 > 
the current is positive irrespective of the value of n. 

Also the magnetic moment is different from the quan- 
tum dot case. The magnetic moment of a particu- 
lar bound state \n,m) is defined as follows: Mn^m = 
q/(2me){n,m\Xz\n,m), where A2 is the z-component of 
the moment of the mechanical momentum A= r x [p + 
(e/c)A]. For convenience we write the magnetic moment 
Mn,m in units of Mq = —eh/{2me), so we obtain 

Mn.m = m + —{n,m\pA^{p)\n,m), (9) 

which is plotted for various one electron states in Fig. ^ 
for the three different systems. In the limit Ba ^ we 
obtain the well known result Mn^rn = 2n -\- \m\ + m + 1: 
the magnetic moment is that of an electron in a homo- 
geneous magnetic field. Notice that all states with the 
same n but to < have the same moment in this limit. 
Furthermore Mn,m > for all bound states with to > 0. 
For Ba ^ 00, the same result as for a circular dot defined 
by hard walls, i.e. Mn^m = is obtained for the models 
with (a) and without (b) overshoot. Thus the magnetic 
moment of the to < states changes sign with increasing 
Ba- This change in sign occurs at a larger i3a-value when 
TO is more negative. Although the extreme limits are the 
same, again the intermediate behavior is different for the 
three systems. As one would expect, the magnetic mo- 
ment in the case of the model with overshoot reaches its 
Ba — > 00 limit at a smaller applied field, than in the case 
without overshoot. For small applied magnetic field we 
also observe oscillations in Mn,m in case of a magnetic 
overshoot. The oscillatory nature smoothes out and dis- 
appears as Ba is raised or when the magnetic field is 
smoother at the edge (see Fig. ^(c)). This behavior can 
be understood by the following picture: depending on 



the radial quantum number n the radial electron den- 
sity has n + 1 maxima. For small magnetic fields the 
electron wavefunction is extended outside the magnetic 
anti-dot. With increasing magnetic field the position of 
the maxima and minima in the electron density shift, and 
when a maximum is at the position of the overshoot, the 
magnetic field has the largest influence, and consequently 
the magnetic moment exhibits a minimum. For a radial 
quantum number n, there will be n maxima in the elec- 
tron density which will shift through the overshoot at 
p/a = 1 with increasing magnetic field strength Ba and 
consequently Mn^m exhibits rt-maxima (one is at Ba=0). 
When a minimum of the electron density is located at the 
overshoot, the magnetic moment has a local maximum. 

In the superconducting disk case the limit Mn,m{Ba 
00) = TO is never reached. In fact, Mn^m slightly increases 
for Ba/ Bq > 35, and we found Mn.m > "th for all values 
of Ba ■ This is a consequence of the small magnetic field 
under the disk, which is always present as one can see in 
Fig. |l|. With increasing applied magnetic field Ba, this 
field grows continuously, which influences the magnetic 
moment. The oscillatory nature of the magnetic moment 
of the superconducting disk case is still vaguely visible in 
the low magnetic field region. 

V. OPTICAL PROPERTIES 

For the present case of magnetically confined dots the 
confinement potential is not quadratic like it is often as- 
sumed for the case of electrically confined dots. This 
has important consequences for the optical absoration 
spectrum. Due to the generalized Kohn's theorerrO the 
long wavelength radiation couples only with the center 
of mass motion of the electrons in the quadratically con- 
fined dots and the absorption spectrum exhibits only two 
peaks. In the present case we observe coupling between 
the center of mass motion and the other degrees of free- 
dom. Transitions are only possible for Ato = ±1 like in 
the case of quadratic confinement, but the other selection 
rule An = 0, 1 is now broken. This is shown in Fig. [l0| 
where the oscillator strength 

f^;^ = {2mjn^){E,,,m - i?o,o)|(0,0|pe±^^|n,TO)|', (10) 

for excitation from the ground state is plotted as function 
of Ba- This becomes more visible in Fig. p] which shows 
an enlargement of the smaller oscillator strength region. 

Consequently the absorption spectrum exhibits a 
larger number of peaks then in the case of a quadratic 
confinement potential. The transition energies as a func- 
tion of the applied magnetic field Ba are plotted in 
Fig. The solid curves correspond to the transitions 
with largest oscillator strength, the dashed and the dot- 
ted curves are respectively for transitions with one order, 
and two orders of magnitude smaller oscillator strength. 
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VI. CONCLUSIONS 

In conclusion we have studied the single particle prop- 
erties of a magnetically confined quantum dot. When this 
is realized through a thin superconducting disk which is 
situated close to the 2DEG of a heterostructure, the mag- 
netic field profile exhibits an overshoot at the edge of the 
disk which leads to a superlinear behavior of some of the 
energy levels as function of the strength of the external 
magnetic field Ba- The consequence of this behavior is 
a discrete filling of the dot as function of Ba- The cir- 
cular current of the electron bound states has opposite 
sign inside and outside the dot region for the m < 
states. These m < states exhibit a transition from a 
diamagnctic to a paramagnetic behavior with increasing 
Ba-&eld. Because of the non quadratic nature of the con- 
finement potential we predict that the optical spectra of 
these new dots exhibits a rich spectrum of lines. 
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FIG. 1. (a) The calculated magnetic field profile in the 
2DEG plotted as function of the radial distance from the cen- 
ter p of the dot for different values of z, the distance of the disk 
above the 2DEG. The inset is a side view of the experimen- 
tal configuration wc have in mind, (b) The vector potential 
profile in the 2DEG plotted as function of the radial distance 
from the center p of the dot for the model systems without 
magnetic overshoot (dashed curve), with overshoot (dotted 
curve) and for the superconducting disk case with z/a = 0.01 
(solid curve) and z/a = 0.2 (long dashed curve). 

FIG. 2. Energy spectra as function of Ba for (a) the model 
without overshoot, (b) the model with overshoot, and c) the 
superconducting disk case with z/a = 0.01. The Sa — > oo 
limiting behavior is indicated at the right of the figure. 

FIG. 3. Energy spectra (in units of ftwc) as function of 
Is/a ~ -s/sT for (a) the model without overshoot, (b) the 
model with overshoot, and (c) the superconducting disk case 
with z/a = 0.01. 

FIG. 4. Radial part R{p) of the wavefunction witli quan- 
tum numbers (n, m) — (0, —1) in case of the model with over- 
shoot as function of the radial distance to the center of the 
dot, for different values of Ba- 

FIG. 5. The filling of the magnetic anti-dot with elec- 
trons as function of Ba for the model with overshoot (solid 
lines) and the superconducting disk model (dotted lines) with 
z/a = 0.01. 
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FIG. 6. The current density profile and the electron radial 
distribution when Ba/By, = 10 are shown for various m- values 

(indicated on the figure) with n = for (a) the model without 
overshoot, (b) with overshoot and (c) for the superconducting 
disk case with z/a = 0.01. 

FIG. 7. The current density profile and the electron radial 
distribution when Ba/Bo = 10 are shown for various n- values 
(indicated on the figure) with m = 1 for (a) the model without 
overshoot, (b) with overshoot and (c) for the superconducting 
disk case with z/a = 0.01. 

FIG. 8. Classical trajectories of electrons confined into a 

magnetic antidot without overshoot for m = — 1 (clockwise), 
m = and m = 1 (counterclockwise). 

FIG. 9. The magnetic moment M„m (in units of 
Af() = —eh/2me) as a function of applied magnetic field 
strength Ba for (a) the model without overshoot, (b) with 
overshoot and (c) for the superconducting disk case with 
z/a = 0.01. 
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FIG. 10. The oscillator strength f^'f^ as a function of the 
applied magnetic field strength Ba for (a) the model without 
overshoot, (b) with overshoot and (c) for the superconducting 
disk case with z/a — 0.01. 

FIG. 11. The same as Fig. |l^ but now the small oscillator 
strength region is enlarged. 

FIG. 12. The transition energy for excitation from the 
ground state (0,0) to (n,m). The solid curves correspond 
to the transitions with highest oscillator strength, the dashed 
and the dotted curves are respectively for transitions with 
amplitude one order and two orders of magnitude smaller. 
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